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HIGHLIGHTS 


• A physics-based model is established to explain capacity fade in Li-ion cells. 

• The model is consistent with the capacity fade trends observed in NCM + LMO/graphite cells. 

• Fatigue crack growth in the graphite may be an important factor in the loss of usable lithium. 
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Capacity fade is reported for 1.5 Ah Li-ion batteries containing a mixture of Li-Ni-Co-Mn oxide 
(NCM) + Li-Mn oxide spinel (LMO) as positive electrode material and a graphite negative electrode. The 
batteries were cycled at a wide range of temperatures (10 °C-46 °C) and discharge currents (0.5C-6.5C). 
The measured capacity losses were fit to a simple physics-based model which calculates lithium in¬ 
ventory loss from two related mechanisms: (1) mechanical degradation at the graphite anode particle 
surface caused by diffusion-induced stresses (DIS) and (2) chemical degradation caused by lithium loss to 
continued growth of the solid-electrolyte interphase (SEI). These two mechanisms are coupled because 
lithium is consumed through SEI formation on newly exposed crack surfaces. The growth of crack surface 
area is modeled as a fatigue phenomenon due to the cyclic stresses generated by repeated lithium 
insertion and de-insertion of graphite particles. This coupled chemical—mechanical degradation model is 
consistent with the observed capacity loss features for the NCM + LMO/graphite cells. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Lithium-ion batteries slowly lose capacity during cycling and 
storage. These capacity losses result from lithium inventory loss 
due to irreversible electrochemical side-reactions such as electro¬ 
lyte reduction, and from loss of the active electrode materials due to 
particle isolation, phase transition, or other degradation mecha¬ 
nism [1,2], Automotive applications require battery systems to 
maintain a sufficient performance for about ten years. Predicting 
the capacity losses that may occur over such a long time frame can 
be challenging. The standard approach is to perform accelerated 
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aging tests where the battery is cycled or stored at a high 
temperature. 

In Part 1 of this work [3] we described the results of an exper¬ 
imental study of capacity fade of NCM + LMO/graphite cells for a 
test matrix covering a broad range of temperature, discharge rate 
and discharge depth. A semi-empirical life model was established 
in Part 1 to account for the measured capacity fade. In the current 
work we build a physics-based model to better interpret the ca¬ 
pacity fade that was measured. This approach can accelerate testing 
by reducing the number and length of battery life tests that are 
needed to extract important data. 

Degradation mechanisms in lithium ion cells are modeled by 
methods ranging from semi-empirical fits of experimental data, to 
complex numerical simulations of the electrochemical side re¬ 
actions. The semi-empirical approach relies upon curve fitting to 
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large data sets and does not necessarily accelerate the testing 
process. In the numerical simulation approach, physical models for 
the degradation mechanisms are combined with the fundamental 
governing equations for the electrode to yield a set of partial dif¬ 
ferential equations. Such simulations require a detailed under¬ 
standing of the specific material chemistry, which can be 
impractical for situations where the chemistry is not well 
understood. 

To create a generalized and widely applicable model, we focus 
upon the dominant degradation mechanisms for composite 
graphite negative electrodes, which are present in nearly all com¬ 
mercial lithium ion batteries. Capacity fade at the graphite negative 
electrode is closely tied to the formation of a solid electrolyte 
interphase (SEI) film on the graphite surface, which can continue to 
grow slowly even after the initial cycles. Current evidence suggests 
that the SEI consists of both a dense, nonporous inner layer on the 
particle surface on the order of 20 nm-50 nm thick, and a porous 
outer layer which has a thickness of up to hundreds of nanometers. 
Growth of the SEI film with extended storage is often modeled as a 
diffusion-limited process, having a thickness growth that is 
approximately proportional to the square root of time (i.e., para¬ 
bolic growth). In most electrochemical modeling studies of calen¬ 
dar loss and cycle loss, it is typically assumed that surface area 
covered by the SEI is uniform and unchanging [4], However, as 
explained next, in situations where mechanical degradation of the 
electrode material is also an issue, changes in SEI area may also 
need to be taken into consideration. 

Electrochemical cycling is known to increase [2] the rate of ca¬ 
pacity loss (i.e., cycle loss) compared to open-circuit or constant 
voltage storage (i.e., calendar loss). Although the volume expansion 
for lithiated graphite is mild compared to lithium alloys, the 
repeated stresses may contribute to cycle loss through mechanical 
degradation at the negative electrode. Intercalation of lithium into 
graphite up to a stoichiometry of LiC6 is accompanied by up to a 10% 
increase in the interlayer spacing. Repeated expansion and 
contraction generates diffusion-induced stresses (D1S), and can 
ultimately result in particle fatigue and crack propagation. The 
stress fields within a spherical electrode particle during lithium 
extraction, as calculated by linear elastic theory [5], are sufficiently 
large to promote fatigue cracking along existing defects. In terms of 
electrochemical capacity loss, useful lithium will be depleted by SEI 
formation on the newly exposed crack surfaces. Particle fracture 
and the electronic isolation of active electrode material are another 
effect of diffusion-induced stresses, although they are not consid¬ 
ered in the present model. 

Recent work has suggested that capacity fade models should 
address both chemical and mechanical degradation [6], For cells 
with a lithium iron phosphate (LFP) cathode, a successful approach 
was to consider mechanical fatigue in graphite electrode particles 
under the repeated DIS during electrochemical cycling [7], Me¬ 
chanical fatigue was then correlated to crack growth on the particle 
surfaces, resulting in Li consumption by SEI formation on newly 
exposed crack surface. Motivation for this model was provided by 
experimental observation that capacity fade in the LFP cells was 
driven primarily by gradual lithium consumption through SEI 
growth on the negative electrode [8], In this present work, we 
follow the approach of Deshpande et al. [7] to better understand 
the effect of both temperature and discharge rate on the capacity 
fade of NCM + LMO cells, guided by insights drawn from a large 
experimental test matrix. 

We present experimental data collected for NCM + LMO com¬ 
mercial power cells containing a mixture of Li-Ni-Mn-Co layered 
oxide and Li—Mn oxide spinel as the positive electrode material. 
The negative electrode contains graphite as the active material. 
Electrochemical cycling was performed in accordance with a test 


matrix which examines the effects of calender life (t), cycle number 
(N), temperature (T), depth of discharge (DoD) and discharge rate 
(/ceil) on the capacity fade. Earlier cycle life and calendar life studies 
for the NCM + LMO cell chemistry covered a considerably smaller 
range of test conditions [9,10] The experimental data for 
NCM + LMO cells presented in Part 1 of this work point to lithium 
loss at the negative electrode as a dominant factor in the capacity 
loss during cycling, similar to results from an earlier study of LFP- 
graphite cells [6], The lithium loss is greatly accelerated at low 
temperature (T = 10°C) and fast discharge (7.5 A, 9.75 A) conditions, 
highlighting the important contribution of DIS to the capacity loss. 
A coupled chemical-mechanical degradation model originally re¬ 
ported in Ref. [7] is modified in the current work to address 
calender loss and variable discharge rate. This allows it to describe 
capacity fade over a wider range of operating conditions. 

2. Experimental methods 

Data was collected on commercially available 1.5 Ah 18650 cy¬ 
lindrical cells produced by Sanyo, a subsidiary of Panasonic. These 
cells have a blend of high capacity Li—Ni—Mn—Co oxide and high 
power Li—Mn oxide spinel as the cathode material, and have 
graphite as the anode material. Extensive details on the test matrix 
and test procedures are provided in Part 1 of this work (Ref. [3]). In 
defining the C-rates, a nominal capacity of 1.5 Ah was assumed. 

Cycling was performed at four temperatures (10 °C, 22 °C, 34 °C 
and 46 °C) and five discharge rates (0.5C, 2C, 3.5C and 5C and 6.5C) 
were employed for cycle life testing. Cells cycled between an SOC of 
1 and 0.5 (i.e., 50% DoD) were selected for modeling in this work. 
The cycle tests were performed between cutoff voltage of 4.2 V and 
2.5 V. During cycling charging was performed at a 2C rate with one 
exception—for the 0.5C discharge samples, charging was per¬ 
formed at 0.5C. During the life cycle tests, cells were stopped at 
weekly intervals for characterization in which the capacity was 
measured by fully charging and discharging the cells at a 0.5C rate. 

3. Model 

3.1. Overview 

Capacity loss at the negative electrode is modeled at the single 
particle level by considering only a graphite particle and omitting 
the other electrode components such as binder, conductive filler 
and other elements of porous electrodes. We assume that lithium 
loss is the capacity-limiting process for the NCM + LMO cells, an 
assertion which is supported by experimental data in the com¬ 
panion paper [3], The depletion of useful lithium occurs through 
the continuous growth of an SEI layer on the graphite particle [11], 
While in reality lithium may also be consumed in repairing the SEI 
layer (which can be damaged by volume expansion of lithiated 
graphite), we have not directly accounted for this in the model. 

Numerical simulations of the SEI layer have yielded a power law 
dependence for the film thickness versus time, L(t) - L(0) = Kt y , 
where y ~ 0.5 [12], The observation that y ~ 0.5 is consistent with 
a diffusion-controlled film growth of the SEI (i.e., parabolic growth). 
Models which take solvent diffusion through the SEI as the rate- 
limiting factor also yield a ~ Vt time dependence for the growth 
of the compact inner SEI layer thickness [13,14], 

In general, test data at elevated temperatures corroborate a Vt 
correlation for calender loss (and presumably SEI thickness) [ 15,16 j. 
Therefore, we model the increase in SEI layer thickness versus time 
by the following expression. 

i-SEl(t) - i-SEl(O) = JCtht 1 / 2 (1) 
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where R t h is a temperature-dependent rate parameter and Lsei(O) is 
the initial SE1 thickness following the cell formation cycle. It is 
emphasized, though, that this expression does not include the ef¬ 
fect of voltage or electrolyte concentration on the SEI formation 
rate. Further, Eq. (1) does not differentiate between SEI formation 
during charge and discharge cycles. One would expect solvent 
reduction on graphite to occur primarily during the charging 
process. 


3.2. Diffusion induced stress 


Symmetric lithium diffusion into, and out of, a spherical particle 
of radius R is described by Fields second law in spherical 
coordinates, 


3C _ D d ( 2 C\ 
at r 2 dr V rj’ 


( 2 ) 


where C is lithium concentration and D is the lithium diffusion 
coefficient in Li x Ce. With an uniform initial concentration of Co, the 
solution for a galvanostatic boundary condition is given by 


c(x, t) = c 0 + 


iR 

FD 


3 2 sin(AnX) 

10 x ^ \^2 sin(%J ) 


(3) 


where x = r/R and r = t/(R 2 /D) [5], The symmetric current density 
over the particle surface is denoted i, where negative currents 
correspond to discharge. The eigenvalues A n are the positive roots of 
tan(A„) = A„. The current density i is estimated from the measured 
cell current by standard expression for sphere packing, 


. I cen R 

3e(AL)’ 


(4) 


where e is the volume fraction of the solid phase within the 
negative electrode, A is the geometric electrode area, and L is the 
electrode thickness. 

For composite graphite electrodes, e ~ 0.65 is a reasonable 
value. Based on the nominal capacity (1.5 Ah) of the NCM + LMO 
cells, and an average energy density of 1.5 mAh cm~ 2 , we calculate 
an electrode area of A = 1000 cm 2 . The negative electrode thickness 
for the cells in this study was measured at 30 pm. 

The evolution of stresses during galvanostatic discharge can be 
calculated from Eq. (3) by the stress-strain relationship for an 
isotropic, linear-elastic solid. The resulting expression in spherical 
coordinates for the tangential stress is. 


1 EQ R 2 I cell 
3 (1 -v) FDa\L 


i(i-rf) 


>y" _g _ 

sin(A n ) 


sin(/ n x) sin(AnX) - (A n x)cos(A n x) 


(5) 


where E is the Young's modulus, v is the Poisson ratio, and Q is the 
partial molar volume of the solute. 

Since full lithium intercalation 6C + Li —► LiC6 results in a 10% 
volume expansion starting from an initial graphite density of 
2.2 g cm -3 , we can estimate an effective average value for the 
lithium solute volume equal to Q = 3.3 cm 3 mol -1 . The Young's 
modulus and Poisson ratio for graphite are estimated as E = 15 GPa 
and v = 0.27. However, the Young's modulus may increase by up to a 
factor of three during Li intercalation [17], A radius of R = 5 pm is 
assumed for the graphite particles in the negative electrode [18], 


A diffusion coefficient of D = 1 x 10 -13 m 2 s _1 is used for Li 
diffusion in graphite between 10 °C and 46 °C. However, Li diffu- 
sivity in Li x C6 is strongly concentration dependent [19], varying by 
almost two orders of magnitude between x = 0 and x = 1. The 
tensile stress, given by Eq. (5), is plotted in Fig. 1 versus dimen¬ 
sionless time at several values of radial position r/R in the graphite 
particle. Negative values for stress denote compression. 

The maximum tensile stress occurs at the particle surface, 
reaching a steady state value given by 


1 EQ R 2 / C eii 
45 (1 - v) FD eAL 


( 6 ) 


after roughly r = 0.2. The cross-over from tensile stress to 
compressive stress occurs at x = 1 /y/2, with compressive stress 
prevailing closer to the particle center. Based on the assumed 
particle radius and Li diffusion constant, the characteristic diffusion 
time constant is R 2 /D = 90 s. In comparison, the time required for 
discharge from SOC = 1 to SOC = 0.5 is 900 s at 2C, and 277 s at 6.5C. 
Therefore, it is safe to assume that the steady state value for tensile 
stress in Eq. (5) is reached each discharge cycle. Values of ae,max are 
plotted in Fig. 2 as a function of cell current / ce ii- 


3.3. Further elaboration on underlying assumptions 

In the current model, graphite particles are treated as isotropic 
spherical particles within linear elastic theory. No distinction is 
made between the basal surfaces and prismatic (edge) surfaces of 
graphite, which have distinct properties in terms of Li intercalation, 
crack growth and SEI layer formation [20], The anisotropic nature of 
graphite is also not included. In choosing our model parameters, we 
try to be consistent with the observed surface morphology and 
porous texture of typical graphite materials used in commercial Li- 
ion cells. Synthetic graphites often have a smooth surface along 
with a moderate mesoporous (2 nm-50 nm) structure, while 
natural graphites have greater surface roughness with little mes- 
oporosity [21 ]. First cycle capacity loss is strongly correlated to the 
graphite surface area, so the N2 BET specific surface areas (SSA) for 
commercial graphites will likely fall between 2 m 2 g -1 and 
10 m 2 g~ J in order to keep first cycle losses below 10% [18], 

We assume that the electrode manufacturing process (e.g., 
calendering) and formation cycle leaves the graphite particles 
with an initial density (p cr ) of sub-micron cracks on surface of the 



as function of dimensionless time, r = Dt/R 2 , in units of MPa. Each trace corresponds to 
the stress evolution at a particular radial position, x = r/R. Positive values of ae indicate 
tensile stress. 
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Fig. 2. Maximum tensile stress at the graphite particle surface plotted versus discharge 
current 


graphite electrode particle [22], To be consistent with the 
observed mesoporosity of the synthetic graphites (i.e, no pores 
<2 nm) the width and length of the surface cracks are presumed 
to have an average value of 15 nm. Similarly we set the crack 
depth to an initial value of 20 nm. We assume that the density of 
cracks on the particle surface is 3 cracks per 1000 nm 2 , which 
works out to pa = 10 6 cracks per 5 pm radius particle. With this 
crack density, approximately 70% of the spherical surface is 
occupied by cracks. 

If we simply consider a graphite particle as a perfectly smooth 
spherical shape, the resulting SSA will be too small (<=d m 2 g -1 ). To 
match values for graphite surface areas, therefore, we choose a 
roughness factor of r gr = 13 for the external surface of the graphite 
particle. The roughness factor in this case is the ratio of actual 
external surface area to the area of a smoothed envelope around the 
particle. By this approach, the initial external graphite surface area, 
excluding the contributions from cracks, works out to approxi¬ 
mately 3.5 m 2 g Cracked surfaces are also assumed to have a 
roughness factor of r cr = 13. Much of this cracked surface area is 
likely to be created during the calendering step of manufacturing 
process or during the formation cycles [22], and therefore would 
not be reflected in the measured SSA's of the pristine material. 
Counting both the external and cracked surfaces of the graphite 
particle, the initial SSA is approximately 10 m 2 g ', which is 
consistent with typical values for fabricated negative electrodes in 
commercial cells. This surface area is distinct from the value used in 
calculating the current density on the graphite particles in the 
negative electrodes, which assumes a smooth spherical surface. 

The chemical composition of the SEI layer is not fully under¬ 
stood and it depends very strongly on the composition of the 
electrolyte [ . It is known to be composite in nature and contains 

both inorganic and organic components. Inorganic compounds 
such as Li 2 C03 are known to be one of the components of SEI films 
formed by reduction of ethylene carbonate. The proportionality 
between SEI thickness and Li consumption is estimated from well- 
known properties of commercial Li-ion cells. In particular, the BET 
N 2 surface area of negative electrodes in commercial cells is 
generally on the order of 10 m 2 g -1 even after electrode 
manufacturing. At the same time, the SEI layer on the negative 
electrode material following formation cycles is in the range of 
20 nm—50 nm [14], and a value of 23 nm is assumed for this 
calculation. 

It is well-known that graphite (365 mAh g 1 ) typically con¬ 
sumes 10% of usable Li on the first cycle. For 1 g of graphite, a Li loss 
of 36.5 mAh (1.36 mmol Li) should be correlated to the formation of 
10 m 2 x 23 nm—0.2 cm 3 of SEI. The proportionality constant is 
given by 


V SE ,= 


1.362 x 10- 3 moles Li 
2 x 10 7 m 3 


= 6000 mol m 3 , 


(7) 


which is considerably lower than the 57,000 mol m -3 for Li 2 C03. 
The explanation for this is that the SEI also contains organic and 
polymeric specifies which do not contain Li. 

Many electrolyte additives such as vinylene carbonate (VC) may 
polymerize to form polymeric components within and over the SEI 
layer [23], Additives are likely used in all commercial cells due to 
their ability to extend cell life. Therefore, the mechanical fatigue 
properties of the SEI film are likely to be different from that of the 
host graphite. Mechanical properties of graphite do not change 
significantly within the narrow temperature range (10 °C—46 °C) 
that is explored in our measurements. On the other hand, some 
polymeric materials have a glass transition temperature within 
that range, and we would expect their mechanical fatigue prop¬ 
erties to change significantly between 10 °C and 46 °C. The me¬ 
chanical fatigue parameters that we include in the model may 
therefore be effective values that reflect an average of both 
graphite and the various SEI components. While mechanical 
degradation of the SEI layer itself is not explicitly included in our 
model, the fitted parameter values may implicitly reflect these 
effects. 


3.4. Mechanical fatigue 

During discharge, the tensile stress on the graphite particle 
reaches a steady-state value given by Eq. (6), and remains there for 
the remainder of the cycle. Lithium insertion during charge, on the 
other hand, produces compressive stresses at the particle surface. 
The repeated fluctuations between tensile and compressive stress 
near the particle surface creates mechanical fatigue, an ideal 
environment for crack growth. Three factors are necessary to cause 
fatigue: (1) a maximum tensile stress of sufficiently high value, (2) a 
large enough fluctuation in the applied stress, and (3) a sufficiently 
high number of cycles. Previous work has demonstrated that 
repeated Li intercalation and de-intercalation is sufficient to induce 
fatigue crack growth in graphite over a long period of continuous 
cycling [5,7 j. 

Paris' Law is a well-known empirical expression for sub-critical 
crack growth caused by mechanical fatigue. 

| = ( 8 ) 

where a is the depth (from topmost surface) of the surface crack, 
and AKi is the difference in mode 1 stress intensity factors for 
charge and discharge DIS. The constants k and m are inherently 
empirical, and can be treated as adjustable parameters in the 
current context. Although we treat them as fitting parameters in 
the model, k and m are related to material properties and should 
be dependent on the temperature, load frequency and other 
environmental factors. An edge crack in a large plate under uni¬ 
axial tensile stress has a stress intensity factor approximately 
equal to 

Ki = (ro,maxb\/na, (9) 

where b = 1.12 for small surface cracks (a <k R). This expression is 
substituted into Eq. (8) to obtain the equation governing crack fa¬ 
tigue growth, 

^ = k(a^ max bV : m) m ( 10 ) 

The solution to Eq. (10) for initial value a(0) — do is given by 
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a = a 0 (^ + 2 ^-k[ae, m ^b^ m afNy m . ( 11 ) 

which describes the crack depth (a) as a function of cycle 
number (N). 

An assumption made in this model is that the number of cracks 
on the graphite particle ( =4nR 2 p a ) remains constant during cycling 
[7], We also assume that the crack length and width remain fixed at 
l cr = 15 nm, as described in Section 3.3. Therefore, only the depth a 
of the surface cracks increases with cycling. The proportionality 
constant between the crack depth and the total surface area 
contributed by the cracks is given by 

A cr = (4uR 2 p cr )(2/ cr r cr ), (12) 

where r cr is the roughness factor. When the crack depth increases 
by an amount da, the newly exposed crack surface A cr da is covered 
immediately by an SEI layer on the subsequent charge half cycle 
with thickness L° EI . Lithium charge consumed per unit volume of 
SEI layer is equal to VseiF where F is Faraday's constant and V SEI 
was defined in Eq. (7) as the molar lithium concentration within 
the SEI. 

The crack depth a increases with each succeeding cycle ac¬ 
cording to Eq. (8), and lithium is lost through SEI formation on the 
newly exposed crack surface. The amount of lithium consumed per 
cycle on the newly exposed crack surface area is then equal to 
\/ SE iFL(! EI (Acrda/dJV). The lithium consumption by the first SEI layer 
on fresh crack surfaces, (dQ/dN)i, is equal to 

(d/v), = K™ xb V^o] m k 

x(l + 2 - m kK ma3£ b^ m ^]v) rai . (13) 


3.5. Calendar loss 

The formation cycle produces an SEI film which we assume to 
cover the entire spherical particle surface with a thickness Lg H . 
With each succeeding cycle, we assume that a new layer is 
deposited on top of the initial SEI film. The total area of the initial 
SEI film is equal to the combined area of the particle surface and 
surface cracks. 

Ai EI =4nR 2 r gr +Acra 0 , (14) 

where r gr is the roughness factor for external graphite surfaces. The 
amount of lithium consumed by the increase in thickness of the 
initial SEI layer during cycle N is given by, 


3.6. Growth of SEI films on the cracked surfaces 

For completeness, it is assumed that the SEI films deposited on 
newly cracked surfaces continue to increase during succeeding 
cycles. (A simpler and less accurate approach would be to assume 
that the SEI grows immediately to its final thickness during the 
cycle in which the crack area is exposed.) If the current cycle is N, 
and we consider a patch of crack surface formed during any pre¬ 
vious cycle i, the SEI layer on that patch will have grown during the 
ensuing N—i cycles. During cycle N, therefore, the total lithium 
consumption by SEI increase on patches formed during previous 
cycles (1 <i<N) is equal to 




(17) 


Referring back to Eq. (13) for (A CT da/dN) and to Eq. (16) for (dLsm/ 
d/V), the expression in Eq. (17) can be written out as 


(§) 3 = W£VkKmaxbv^] m 

x (l -I—2— k[(Tg tm3X bV^r\ a Q 2 i^ 

x^(iv-r 1/2 (^). (18) 


This term reflects the coupled effect of chemical and mechanical 
degradation. 

The total capacity loss per cycle for a single graphite particle 
equals the sum of Eqs. (13), (16) and (18), 



However, for comparison to the measured capacity loss, Eq. (19) 
must be expressed as a fraction of the nominal capacity of the as- 
received commercial cells. To calculate the fractional capacity loss 
predicted by the model, the initial capacity of a single graphite 
particle (following formation cycle) needs to be estimated. 
Assuming a formation cycle efficiency of 7jf C = 0.9 and a graphite 
specific capacity of 365 mAh g _1 , the capacity at the start of cycling 
is 


Qo = V fc (365 mAh g- 1 x |uR 3 Pgr ), (20) 

where p gr = 2.25 g cm 3 is the density of graphite. 


4. Results and discussion 


Qsei(N) = V SEI FA° EI [I SEI (N) - F SE .(iV - 1)], (15) 

where dlsEi/dN ~ Lsei(FI) - Lsei(N - 1) is the incremental change in 
SEI thickness for cycle N. Based on Eq. (1), the increase in SEI film 
thickness follows a square root of time dependence. The amount of 
lithium consumed with each successive cycle is therefore 

Eq. (16) describes the total lithium loss per cycle, (dQ/dN)2, due 
to increase in thickness of the original SEI layer. The time varies 
approximately linearly with cycle number, and (dt/dN) is therefore 
the slope of the elapsed time versus total cycles. 


The cumulative capacity loss versus cycle number, as predicted 
by the model, was determined by numerically integrating Eq. (19) 
(which gives the capacity loss per cycle). The Paris' law constants 
k and m, and the SEI layer rate constant R t h were treated as 
adjustable parameters, and were fitted to measured data by a 
constrained non-linear optimization routine in MATLAB. A single 
global value of m = 2.2 was used for all temperatures and discharge 
rates. On the other hand, k and Kth were allowed to vary with 
temperature. The fit results for the 10 °C and 46 °C data are shown 
in Figs. 3 and 4, respectively. Fitted values for m, k and K t h are 
summarized in Table 1. The experimental NCM + LMO capacity loss 
is normalized with respect to the cell capacities measured before 
the first cycle (Qo). 
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Fig. 3. Capacity fade data for NCM + LMO cells cycled at 10 °C and discharge rates between 0.5C and 6.5C. Solid black lines indicate the total capacity loss calculated from Eq. (19), 
while the blue, red and green lines represent the relative contributions of the crack growth (Eq. (13)), SEI layer growth (Eq. (16)) and SEI layer growth on cracked surfaces formed 
during previous cycles (Eq. (18)), respectively. The model is calculated from parameters listed in Table 1. (For interpretation of the references to color in this figure legend, the reader 
is referred to the web version of this article.) 


As shown in Fig. 3, the dominant lithium loss mechanism at 
low temperature and high discharge rate is crack propagation, 
(dQ/dN)i, combined with SEI formation on newly exposed 
graphite surfaces. At 46 °C, however, growth of the existing SEI 
layer, (dQ/dN) 2 , appears to be the dominant mechanism for ca¬ 
pacity loss—consistent with commonly observed behavior. A 
closer examination of the SEI thicknesses and surface crack 
depths predicted by this model provides a coherent picture of 
negative electrode degradation during extended electrochemical 
cycling. 


Based on the assumption of 90% formation cycle efficiency, the 
initial SEI layer thickness is Z.° Ef = 23 nm. The final SEI thickness 
varies depending on the temperature and discharge current. For 
cells cycled at 46 °C/0.5C conditions, the predicted SEI thickness is 
100 nm after roughly 1900 cycles. On the other hand, at 10 °C/0.5C 
cycling the predicted SEI thickness is only 60 nm after more than 
3500 cycles. In comparison, experimental measurements of the 
initial SEI thickness in LMO/graphite cells were on the order of 
40 nm, which increased to nearly 400 nm after storage for 1 year at 
40 °C [14], 



while the blue, red and green lines represent the relative contributions of crack growth (Eq. (13)), the SEI layer growth (Eq. (16)) and SEI layer growth on cracked surfaces formed 
during previous cycles (Eq. (18)), respectively. The model is calculated from parameters listed in Table 1. (For interpretation of the references to color in this figure legend, the reader 
is referred to the web version of this article.) 
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Table 1 

Fitted values for model parameters. 

T°C fc/10~ 20 Kth/10 -10 m 

10 13.6 16.2 2.2 

22 3.9 18.2 2.2 

34 2.5 25.4 2.2 

46 3.6 45.1 2.2 


In Section 3.3 we assumed that there was a fixed population of 
cracks (10 6 ) on a single graphite particle (possibly formed during 
manufacturing and formation cycling), each having a uniform 
initial crack depth of 20 nm. The increase in crack depth that is 
predicted by the model, as expected, is largest at 10 °C and at high 
discharge rates. At 10°C/6.5C conditions, the final crack depth is 
65 nm after 1000 cycles. In contrast, the crack depth does not 
change appreciably up to 2000 cycles for cells cycled at 0.5C rate at 
any of the measured temperatures. 

Based on the values of initial crack depth (ao), particle size (R) 
and initial crack density (p cr ) listed in Table 2, the initial 
graphite surface area is roughly 10 m 2 g -1 (see Section 3.3). For the 
10°C/6.5C cycling conditions (which generated the deepest cracks) 
the surface area is more than doubled, with a value of 25 m 2 g -1 at 
1000 cycles. These surface areas are within the range observed for 
synthetic graphites [20,18,2- ;. 

Of the three adjustable parameters included in this model, two 
of them (k, K t u) exhibited variation with temperature. Fitted values 
of the fatigue parameter k are plotted versus temperature in Fig. 5. 
The value of k is on the order of 1 x 10 20 , increasing sharply as the 
temperature is lowered to 10 °C. This strong temperature depen¬ 
dence is inconsistent with the known mechanical properties of 
graphite, which does not show any transitions between brittle and 
ductile phases near 10 °C. 

One possible explanation for this anomalous temperature vari¬ 
ation of k is that it is implicitly parameterizing other phenomena, 
such as the mechanical degradation of the SEI layer itself. If the SEI 
layer were to deform and crack during electrochemical cycling, 
usable Li would be consumed in repairing the SEI layer. As noted in 
Section 3.3, many electrolyte additives such as VC can be reduced to 
form polymer-like components embedded within the composite 
SEI structure. Since certain polymers have a brittle-to-elastic glass 
transition temperatures around 10 °C, and consequently exhibit 
large changes in the mechanical fatigue behavior, this would be 
consistent with the steep rise in k at that temperature. Based on the 
information we currently have, it does not appear possible to 
separate the cracking of graphite active material from cracking of 
the SEI layer. 


Table 2 

List of parameters. 


Symbol Description Value 


A Geometric area of electrode 1000 

Do Initial crack depth 20 

D Li diffusion coefficient 1 x 10 -13 

E Young's modulus of graphite 15 

L Electrode thickness 30 

l a Crack length, width 15 

m Paris' law exponent 2.2 

R Radius 5 

Vu Li concentration in SEI 6000 

e Solid phase porosity 0.65 

i] Formation cycle efficiency 0.9 

p a Cracks per 1000 nm 2 3.18 

Q Partial molar volume of solute 3.3 

v Poison ratio 0.27 

Graphite capacity 365 


Units 

cm 2 

m^s - * 1 

GPa 

pm 


mol m 3 



Temperature / °C 


Fig. 5. Fitted values for the Paris' law constant k plotted versus temperature. Solid line 
is included as a guide. The +x error bars indicate the maximum temperature excur¬ 
sions of the cells during 6.5C discharge. 


The rate constant for SEI layer growth, K t h, is plotted versus 
temperature in Fig. 6. The temperature dependence of K t h is 
consistent with the Arrhenius equation, with an activation energy 
£ a = 21.2 kj mol -1 and pre-factor K'th.o — 1.1 x 10 -5 . This is consis¬ 
tent with a thermally activated chemical reaction rate for SEI for¬ 
mation. Further, this activation energy is similar to the value 
(24.5 kj mol -1 ) that is obtained in Part 1 of this work by using an 
empirical fitting approach [3]. 

A major assumption made in this analysis is that the cell tem¬ 
perature remains equal to the temperature of the environmental 
test chamber. During rapid discharges, however, the internal tem¬ 
perature of the cells can rise by 20 °C or more, followed by a cool¬ 
down when the battery is recharged at a moderate rate. To test this 
assumption we monitored the external battery temperature during 
6C discharge. The results (shown in the Supporting information, 
Fig. SI1) indicate that cells in a 22 °C environment can rapidly 
heat to 50 °C during a 6C discharge. Similarly, cells in a 50 °C 
environmental chamber heat up to 75 °C during 6C discharge. The 
internal cell temperatures should likely exhibit even larger 


t r c 



temperature. The solid line represents the best-fit to the Arrhenius equation. The -x 
error bars indicate the maximum temperature excursions of the cells during 6.5C 
discharge. 
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temperature spikes than those measured on the outer surface by 
thermocouple probes. 

The primary impact of these cyclic thermal fluctuations should 
be on the reaction kinetics of the SEI formation, which would speed 
up during the high T intervals. Therefore, the f<th values listed in 
Table 1 may actually correspond to effective temperature-averaged 
values. In Figs. 5 and 6 we include error bars in the +x direction to 
indicate the maximum temperature excursions experienced by the 
NCM + LMO cells (at 6.5C). 

Agreement between NCM + LMO cycling data and the chem¬ 
ical-mechanical degradation model provides strong, if indirect, 
support for the importance of fatigue cracking and mechanical 
degradation towards overall cell capacity fade. The validity of this 
model hinges upon gradual changes in the surface morphology and 
crack growth of graphite particles in response to extensive elec¬ 
trochemical cycling. To date, however, there do not appear to have 
been any quantitative studies of graphite cracking and surface 
morphology during repeated Li intercalation. 

In Part 1 of this work [3], we provide a detailed differential 
voltage analysis that supports our assumption that lithium loss is 
the dominant cause of capacity fade in the NCM + LMO cells. From 
the shifts and spacing between the characteristic dV/dQ peaks of 
lithiated graphite, it is determined that Li loss outpaces both 
cathode material loss and anode material loss as the capacity- 
limiting factor in almost all the tested NCM + LMO/graphite cells. 
However, this provides only indirect support for fatigue cracking 
and coupled chemical-mechanical degradation. Direct imaging 
studies of graphite particles [25] and cathode materials [26] 
following electrochemical cycling have found evidence of crack 
formation and propagation, but there is no way to correlate this to a 
change in the overall graphite surface area. At the same time, 
previous results have demonstrated that the BET surface area of 
graphite particles can be increased more than four-fold by me¬ 
chanical milling [20], In this context, the two-fold increase in sur¬ 
face area through fatigue crack growth predicted here appears 
plausible. 

5. Conclusion 

Capacity fade in lithium-ion battery with a blended 
Li—Ni-Co-Mn oxide and Li-Mn spinel positive electrode was 
studied in cells cycled across a range of temperatures (10 °C-46 °C) 
and discharge currents (0.5C—6.5C). The experimental data was 
described by a single particle model which combines two mecha¬ 
nisms for lithium loss in the NCM + LMO cells: (1) fatigue crack 
growth at the graphite particle surface from repeated diffusion- 
induced stresses, (2) growth of the SEI film thickness under 


diffusion-limited chemical kinetics. This dual degradation model 
accurately predicts capacity loss across wide ranges of operating 
temperature and discharge current, while simultaneously 
providing an intuitive physical picture of degradation at the nega¬ 
tive electrode. 

Appendix A. Supplementary data 

Supplementary data related to this article can be found at http:// 
dx.doi.org/10.1016/j.jpowsour.2014.07.028. 
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